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Abstract 

We prove a set of identities for the anomalous dimensions of the quark and gluon conformal 
operators in the flavour singlet channel in QCD. These relations arise from the graded commutator 
algebra of the M = 1 superconformal group. We evaluate the rotation matrices for the quantities 
under study from the conventional dimensional regularization to the supersymmetry preserving 
regularization scheme. Using them we verify the equalities in two-loop approximation employ- 
ing the results for the NLO anomalous dimensions of the conformal operators in the minimal 
subtraction scheme derived earlier. 
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1. Introduction. In spite of a number of attractive theoretical features of supersymmetric gauge 
theories |l| in view of unification of the fundamental forces of Nature, the manifestation of their 
predictions was not observed experimentally so far. Nevertheless, they provide an excellent tech- 
nical playground for the exploration of new physical concepts. It is remarkable that due to the 
high symmetry of the underlying Lagrangian the theory enjoys the property of reducibility of 
independent parameters such as gauge couplings and field renormalization constants^. Since the 
structure of the corresponding Lagrangians resembles that of ordinary gauge theories with an ad- 
justed fermion sector it can serve as a technical tool for deriving relationships between observables. 
Using heuristic arguments of this kind several relations have been derived in Ref. between the 
anomalous dimensions of the local quark and gluon operators without total derivatives which ap- 
pear in the description of the deep inelastic scattering via the operator product expansion. One 
of them is the empirically established Dokshitzer relation ||. 

In this paper we address the issue of constraints on the anomalous dimensions of the conformal 
operators from the point of view of graded commutator algebra of the superconformal group and 
Ward identities for the Green functions with composite operator insertion. Presently we consider 
restricted Q-supersymmetry transformations which provide relations for the scale anomalies of 
composite operators with total derivatives, i.e. their anomalous dimensions. We make use of 
the subalgebra of the full superconformal algebra ||, which includes scale transformation and 
consists of the relations, apart from the usual ones, 

[Q,V,]- = 0, [Q,M, u ]- = ^a, u Q, [Q,Q] + = 2^V„ [Q,V]_= t -Q. (1) 

In the infinitesimal form their action on the field operator <p is defined as 5 G <p = i[4>,Q}- with 
Q = Vfj,, M.fj.vj'D for even generators and Q = (Q for odd Q with ( being a Majorana Grassman 
valued spinor. One can equally include special supertransformations [f|, |5| which would allow to 
derive relations for the special conformal anomalies by means of the anomalous superconformal 
Ward identities. 

2. M — 1 QCD. The Lagrangian for the M = 1 super- Yang-Mills theory M in the Wess-Zumino 
gauge H takes the form which resembles ordinary one-flavour QCD with Majorana fermions, 
ip = Cip T , in the adjoint representation of the SU(N C ) group 

£d = -\ (G« u ) 2 + % -r + \ {D a f • (2) 

The Wess-Zumino gauge breaks linear M = 1 supersymmetry but the remaining short supermul- 
tiplet ip a , D a ^j respects the restricted non-linear supersymmetric transformation laws 

$ V = "G^cvC - ^ a 7sC, 5 Q Bl = -iCi,r, 5 Q D a = CT b l^ b - (3) 

2 The latter holds true provided a supersymmetry preserving quantization and regularization procedures are 
used to handle the underlying Lagrangian. 
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The action transforms w.r.t. these transformations as 5®S = — J d 4 x ^d p C,Q p — X)</> (SS/8(j)) > 
with the non-anomalous (on quantum level) || supersymmetry current Q p = \G a ^ u o ^7 p ip a f pT0|1 . 
It is conserved on the mass shell: d p (Q p = J2<j, (^S/8<p) 6 Q (p and, thus, 5 Q S = 0. 

Eqs. (^|) form, however, a modified algebra which contains apart from the usual terms on the 
r.h.s. of ([I]) also a gauge transformation #s au § e j e = —210^5^ + <5 gauge , with field 

dependent gauge parameter, 2i(i B a ( 2 , and translation vector a p = (11^(2- Note, however, that 
they define the ordinary SUSY operator algebra (|l|) in the space spanned on gauge invariant 
objects. 

In our consequent discussion we will be most interested in the commutator algebra restricted to 
the action on the "good" light-cone components [[EJ of the field operators introduced as follows^ 
with projectors IJ± = ^7^7± for fermion fields: the i-components of which are defined via 
if)± = II±ifj a ; and with the two-dimensional metric tensor g^ v = g^ — n^n* — n u n*, for bosons: 
B a ^~ = g-^ u B®. We also impose a restriction on the constant Majorana fermion £ + = i7+C = 0, 
so that 5®B+ = 0. Then a peculiar feature of the supersymmetry transformations (|3]) becomes 
manifest, namely, that they do not involve the auxiliary field, D a , provided we restrict ourselves 
to the "good" components — the ones which enter the quasi-partonic conformal operators defined 
below in Eq. (§). Explicitly |3|, |4 



S Q r + = -Gfry-T^C, ^Bf = -Klin- (4) 

Moreover, they form a representation of the unaltered supersymmetry algebra ([!]). 
3. Supermultiplet of conformal operators. To derive the relations for the anomalous dimensions we 
have first to find an irreducible representations of the superalgebra in the basis of the conformal 
operators 0: 

(5) 

where d = d+d and T>=V — T>, and the factor ^ in front of the fermion operator serves to avoid 
double counting of the same components of the spinors due to their Majorana nature. 

According to the appendix one can introduce the following combinations of the conformal 
operators 





j 3l Jl> |p2[ j _|_ 1 01 J + l J" 



3 The + and — components of any vector are obtained by contraction with the two light-like vectors n and n* , 
such that n 2 = n* 2 = and nn* = 1. 
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with r = V(A) for the upper (lower) entry on the l.h.s. of these equations. By construction they 
transform covariantly w.r.t. supertransformations and together with 

(j+2)(j + 3) 

U + i) + " 




-Gi^d + ypr\^\it{ (7) 

ji 



form an irreducible representation with the transformation laws: 

s Q sj t = \[i - (-iy\ (Vj-u, 5 Q 4 = ~[i - (-i)i cv jh (8) 

S Q Vj l = ~[1 + (-iy\ Qtj-u, = ~[1 + (-1)'] CM,-,, (9) 

^Vi-u-j = -T-C {Sj t + S*_ lt } - 7-7sC + P?_ u } • (10) 

The transformation law for the operator Uji follows from the observation Uji = — 75 Vj; and the 
Eq. (0). 

4. Ward identities and commutator constraints. To proceed with our derivation of the Ward iden- 
tities we have to fix the remaining ordinary gauge degrees of freedom of the classical Lagrangian 
(HD, i.e. to add a gauge fixing term together with an associated ghost piece. Provided we would 
work in the superfield formalism this can be done in a way which preserves linear supersymmetry 
with SUSY Fermi-Feynman gauge. However, then we should proceed with the complete gauge 
supermultiplet at an expense of a number of auxiliary fields on top of the dynamical ip and B. 
To achieve this goal in the Wess-Zumino supergauge without explicit breaking of the supersym- 



metry on the Lagrangian level one is forced to use the light-cone gaugen [fT3|, |19|, 14 1 so that 
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C = £ c i + £ g f + £ g h, where £ gf = — ^ 1 £gh = oo a V+uj b p0[ . Ghosts decouple from all Green 

functions since n fl {B^(xi)B^(x2))\^- f o — 0. However, because the Lorentz invariance is violated 
by fixing a particular direction in the Minkowski space with the vector n M , different (i.e. "good" 
and "bad" ) components of the field operators renormalize with different renormalization constants 
and, moreover, non-local counterterms are required to cure all divergencies of the classical action, 
i.e. [EDI 



b; - z l J 2 (b; - (1 - z, l )n,n a ) , v - z\ /2 (n + + z 2 n^) v>, g - z gg , z - z 3 ^, (11) 



4 Covariant gauges inevitably break SUSY since the necessary condition for keeping symmetry of the Lagrangian 
intact is that the algebra of the gauge and rigid symmetries has the form of a semi-direct product [fl5|| . However, 
this breaking is only due to the BRST exact operator S BRST 6 Q (oj a d f _ l B^) (since [JBRST^ 5 Qj_ = @ @) and 
will not affect physical quantities. We shall address this question within the present context elsewhere. Next, the 
axial gauge breaks the supersymmetry of the Lagrangian but leads to supersymmetric counterterms for the gauge 
multiplct However, the latter property is violated as soon as matter superfields are taken into account |is}| . 

In spite of the fact that we are considering only gauge multiplet we prefer to deal with supersymmetry preserving 
gauge fixing on the Lagrangian level. 
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where Q a = \T>+ 1 G + -j is the only non-local structure involved. This instance complicates the 
derivation and the use of the conformal Ward identities. However, the fact that only "good" 
components are relevant suggests to integrate out the "bad" ones in favour of the "good" ones 



13, 14, in the path integral for the Green functions with operator insertions, O, 



(OX) = lim (J D(j) e iS j J D<j) OXe iS , (12) 

where S = J d d xC and X is a monomial of the "good" components of the field operators, i.e. 
X = Hi 4>+( x i) rij i>+( x j) Ilfc B ± (xk). In spite of the fact that this procedure results in the non- 
local form of the action, the Lagrangian manifests explicit supersymmetry w.r.t. renormalized 



transformations (|4]) since from the Slavnov- Taylor identities it follows that Z 3 = Z 2 = Z$ [[14 
To maintain this property in perturbation theory it is necessary to deal with UV divergencies 
in a supersymmetric way, i.e. by regularization via, e.g. Siegel's dimensional reduction^] (DRED) 
|2"2"| , |2"B| , |2~4| . Provided we have done this, the light-cone gauge preserves SUSY (ffi) on the quantum 
level as well and the transformations @ remain unaffected even for the renormalized fields. This 
was checked in one-loop approximation for the principal value (PV) |19| and the Mandelstam- 
Leibbrandt (ML) [|14[] prescription on auxiliary ^q— pole in the gluon propagator. In the latter case 

— 1/2 

Z g = Z 3 , while in the former Z 2j 3-factors are e-dependent due to the fact that the PV prescrip- 
tion violates power counting. This can be traced back to the breaking of rescaling invariance of 
the gluon density matrix. 

It is obvious, that in the course of renormalization the operators S" 1 (V 1 ) mix with each other. 
We introduce the renormalized operators according to 

lO^± {Zohk Z^, with Z =(J5), *,= (f z .?), d3) 

where Oji stands for a two-dimensional vector «S = (52), an d similar for V % . The superscript (0) 
means that the operator is written in terms of bare fields. And the anomalous dimensions are 
defined via the renormalization group equations 

^[Oii] = -E7S[0«], with 7°=U 7 o^o| for = S,V, (14) 



k=0 



^N = -EVN. (15) 

a h L k=0 

The fermionic operators V and U evolve with the same anomalous dimensions, A^, since Uji = 

— 75 Vj/. Note, that from symmetry properties of the operators it follows that the anomalous 

5 This is not a fully self-consistent regularization method |24j but at least its reliability as a scheme which respects 
supersymmetry has been checked on two- loop level [E3| . 
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dimensions for the parity even/odd operators vanish for even/odd j, while they admit both values 
for fermionic case. 

For the derivation of the Ward identities we proceed with the regularization of the action via 
dimensional reduction since otherwise 5®£ c i = —^f abc ip a 'j^ip b C'j fl ip c due to inapplicability of the 
4-dimensional Fiertz identity. Although it is a legitimate procedure to use a symmetry breaking 
regulator, but this would result to an addendum in the Ward identity and as a result in the 
relations we are going to derive. With Siegel's method at hand the supersymmetry Ward identity 
simply reads 

([o^x) = -<m-^>, (is) 

with S Q Oji given by Eqs. (|^)- ([T0|) . Finally, the dilatation Ward identity |J adjusted to the present 
scheme is 

{[On\8 D X) = E{(/ + 3)1 + 7 °L ([Ou)X) +J2^(g)(i[O j A i ]X), (17) 

k=0 i 

where the last term stands for Green functions with the renormalized operator insertions, i.e. 
differential vertex operator insertion [A g ] = g-§^£ and equation of motion operators i?^ = (j>j^- 
Here 1 + 3 = 1 — 1 + 2{d C Q n + 1) = I + 2d < Q n , for gluonic and fermionic operators, respectively, with 
canonical dimension^] c?^ an of the field 0. Thus from the last commutator^ in Eq. (|I|) acting on 
the Green function with conformal operator insertion (12) ([Oj;][<5^, 8 D ]_X) = ^([Oji\8®X) and 
from the Ward equalities displayed above we can easily obtain a set of identities for the anomalous 
dimensions. For this one has to note that ([OjiA g ]X) = g^([Oji]X) and J2(f>=tji,G l<t>{[^ji^A^) = 
i^^N {[O ji]X) with 7^, = 7g = l<f> due to the gauge we have chosen. N is the total number of 
fields in X. Therefore, these terms drop out from the commutator of the scale and supersymmetry 
transformations and thus do not show up in the constraints. Namely, for O = S we have in 
components 

n n n—1 

E {^n+l^m+lD^m,/] + 12 72n+l,2m+l D^m+l,;] } = E ^2n,2m [V 2m ,i] + E ^2n,2m+l [V 2m +l,z] , 



m=0 m=0 m=0 

^ 2m+l,/J , 



E { 2 ^y2n+l,2m+l[^2m,l) + 22 72n+l,2m+l D^m+l,/] } — E ^2n+l,2m [V 2m j] + E ^2n+l,2m+l [V; 

m=0 m=0 m=0 

while for O = "P 

n—1 n n—1 n—1 

E n 72n, 2m+ 2[W2 m+ l,/J + E 12 72n, 2 JW 2m ,,,] = E A2n-l,2m A2n-l,2m+l[^2m+l,z] > 

m=0 m=0 m=0 m=0 

n—1 n n n—1 

E 72n,2m+2[^2m+l,«] + E 72n,2m 2n,2m ^2n,2m+l[^2m+l,i] • 

m=0 m=0 m=0 m=0 



6 Obviously, only "good" components are endowed with well denned canonical dimensions. 
7 Note that its nonzero r.h.s. just serves to shift the canonical scale dimension of the elementary field <j> by \ in 
mass units, thus changing, for instance G?£f n — > dg n + 5 = (i™. 
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Then, using linear independence of conformal operators we get finally the following relations 



V 

/2n+l,2m+l 
/2n+l,2m+l 

/2n+l,2m+l 

22 $ 
7 2 n+l,2m+l 

72n+l,2n+l 



/2ra,2m 
/2n,2m+l ~~ 

v 

/2n+2,2m ~~ 

V 

/2n+2,2m+2 



A 2 n,2m, m <n, 



^2n,2m+l, 
^2n+l 

= A 



to < n — 1, 
m < n , 

2n+l,2m+l; VTL < 11 







V - 

72n,2n — u ■ 



(18) 
(19) 
(20) 
(21) 
(22) 



This is the most general set of equations for the anomalous dimensions of operators with total 
derivatives which arise in a theory with M = 1 supersymmetry. 

5. Reduced supersymmetry relations. From the above relations for the anomalous dimensions of 
the operators which form a supermultiplet we can deduce identities for the familiar anomalous 
dimensions of quark and gluon operators. For this purpose note that they are related to each 
other via the following matrix equation 



( l+,*\ 

fc+1 hk 
k 'Jk 
V fc+1 "^jk J 



( 



2fc + 3 



V 



1 


fc+3 


6 


fc+3 


\ 




6 


j 


J 


-1 


A: 


6 


fc 




6 


3 


3 






(fc+3)(j+3) 


6 


fc+3 




j+i 


60+1) 


3+1 


i+i 




i+3 


k(j+3) 


6 


fc 


/ 


3+1 


6(i+l) 


J+1 


i+l 



QG, 



Ijk 

V GG 7,fc J 



(23) 



Then the transformation to the usual operator basis becomes trivial. Let us consider some par- 
ticular limits of Eqs. (]18|)-(|2~T|). Namely, from Eqs. (|22| ) we obtain the well-known Dokshitzer 
relations empirically established in the original paper || for the vector channel (here and below 

In = 7i) 



QC 



+ 



6 



G % = { Q % + G % 



J 6 

From (1181) and (Elf), reduced to the forward case, we have 



Q( 



y ^ + j+T 



7 j+i 



QQ, 



7 ' ~ 6 



J -Q G jf and QQ jf +1 + 



V,A. 



1 3+1 



(24) 



(25) 



j + 1 " J 6 

Finally, from Eqs. ( |T9"D and (^0]) a relation follows between the diagonal (read forward) and non 
diagonal elements of the anomalous dimensions of the conformal operators 



" i 

J 



1 + 3 QGV 



A 



6 

j 



7,- 



J + 1 
2j + l 



^j+i,j-i> 



' i 

J 



■7+3 
6 



QG, 



7, 



J + 1 



2J + 1 



A i+i,i-i> (26) 



i+ij-i 



GG' + ^ - QC V 

'J+l,i-l T a 1 3+1,3-1 



6 



1 3+hJ-l 



QQ, 



Ij+lj-V 



(27) 



j + 1 6 j + 1 

Obviously, in the leading order conformal operators do not mix and the r.h.s. of the equations 
(f26|) are zero. However, beyond one loop || these equations provide a non-trivial check of existing 
results. 
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Figure 1: Generic form of Feynman diagrams for z projected onto the tensor structures given in 
the text. 



6. Transformation from SUSY preserving to conventional DREG schemes. Let us address now 
the question of explicit checks of the above predictions beyond leading order of QCD perturbation 
theory. Since all above equations hold only for the entities evaluated by means of supersymme- 
try preserving dimensional reduction we have to compute rotation matrices to the conventional 
dimensional regularizationP] (DREG) — a scheme used in practical QCD calculations where all 
higher order results are available. Note that the supersymmetric limit of ordinary QCD can be 
achieved by equating the Casimir operators Ca = Cp = 2Tp = N c . 

The change of the scheme is achieved via the following finite transformation of the quark- 
gluon operator O = Cgq) renormalized according to the conventional DREG to the DRED 
scheme [0] DRED = z[0] DREG . Thus, the quark-gluon anomalous dimension matrix, 7, for the 
regularization with DRED are related to the DREG one via 

7 DRED = Z7 DREG 2 -1 _ f3( g )JL z . z ~\ ( 28 ) 

dg 

In order to check the relations we have derived above in the two-loop approximation the 
problem is thus reduced to the computation of z at 0(a s ) (see Fig. [1]). Let us add few remarks 
on this calculation which has been performed with ordinary QCD Feynman rules identifying 
afterwards the Casimir operators. Since the Clifford algebra is considered as 4-dimensional in 
the DRED as well as in conventional DREG the projectors used are the same in both schemes 
Q^iV'A) _ k^j 5 ^j_ f or the parity even and odd sectors, respectively. On the other hand the 
gluon polarization vectors are treated as 4-dimensional in DRED and <i-dimensional in DREG. 
We have then for parity even case G 7r — (d — 2)~ 1 <r^ d ) with d — 4 for DRED and d = 4 — 2e 
for conventional DREG, while for odd parity g tt^ u = i(d — 2)~ 1 (d — 3)^6^-+- We have used the 
HVBM scheme |32] for dealing with 75 and e^ uprj which are pure 4-dimensional objects. These 
procedure has proved to be the most reliable for these purpose. Moreover, as a cross check we 



have adopted also Larin's prescription |33| according to which we have used the substitution 
7aj75 — — ^tivpalvlpla m quark-gluon diagrams and the resulting product of two e-tensors has 
8 Remarkable that the number of different transformation matrices in the forward case equals to the number of 
independent investigations on the subject Diversity of opinions is welcome but outside physics. 
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been understood as e IMUp(T e al3 "' s = —Algj^gPgJgf) with the metric tensors being <i-dimensional for 
DREG and 4- dimensional for DRED. Note, however, that we have been able to obtain the same 
results as for the HVBM recipe only for QG and GG graphs and have failed for others. Due to 
gauge invariance of the rotation matrices we have done the calculations using covariant Feynman 
and non-covariant light-cone gauges with indeed identical final results. 

Finally, we have found the finite part of the anomalous dimensions of the quark and gluon 
conformal operators, obtained from the difference of the supersymmetry preserving and the con- 
ventional dimensional regularization schemes, which reads 

z jk = 16 jh + ^N c {zf5 jk + + (-l) j - h )} , (29) 

with the following diagonal matrices 

/ j(j+3) 12 \ / j(j+3) 12 \ 

Z V,V = 2(j+l)(i+2) j(j+2)(i+3) ^ Z D,A _ 2 (j+l)(j+2) j(j+l)(j+2) ^ (30) 

V 6(J+2) _ 6 / V ~ 3(j+l)(jf+2) ~6 _ (j+l)(i+2) / 

for the vector and axial channels, respectively, and the universal non-diagonal part 







6(2A;+3) 



ND " fc(fc+l)(fc+2)(fc+3) I /qi\ 

i k I (2fc+3) (2fc+3)(i-fc)(j+fc+3) I ' X ' ■ ' 

V 6(fc+l)(fc+2) fe(fc+l)(fc+2)(fc+3) 

In the course of the calculation we have clarified the reason for a difference in the rotation 
matrices given in the literature for forward scattering. For instance, the recent results of Ref. 



[pl| can be reproduced in practically all cases (except axial GQ channel) provided we discard 
the contributions of e-scalars in the external lines. We can hardly advocate this recipe since the 
z-matrices have to be considered as insertions into the internal virtual lines, so that the e-scalars 
contribute on equal footing with the <i-dimensional gauge particles. The diagonal eigenvalues, z®, 
coincide with |2^, [28| for parity even and with [|30| for parity odd operators. 

Note that only diagonal elements of these matrices are required to fulfill the supersymmetry 
relations in Eqs. (^4|) and ( p5|) which can be explicitly checked transforming the results of Refs. 
[p9 , j30| , j3~4j] with the help of the rotation matrices Cj30|) to the DRED scheme. In Eq. ( |2~6|) the off- 



diagonal elements 0] enter which have to be rotated with (|3l|) to the supersymmetry preserving 
scheme as well. The appearance of these turns out to be the reason why the authors of Ref. j31] 
have found the violation of the third SUSY relation in NLO — they have used an equation of [Q] 
which corresponds to our Eq. fl2HD but with zero r.h.s. 

Moreover, we have found that the general relations ( p^ - |2l|) are satisfied when transformed with 
the rotation matrices we have derived presently. This provides a confirmation for the correctness 
on the two-loop anomalous dimensions for the QCD composite operators available in the literature 
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7. Summary. To summarize we have derived in the present study a general set of relations for the 
anomalous dimensions of the bosonic and fermionic conformal composite operators in the Yang- 
Mills theory with M = 1 supersymmetry. In their reduced form two of them were known earlier 
while the last one (Eq. (p6|)) contains a novel non-vanishing r.h.s. with the non-diagonal elements. 
This explains the difficulties in the check of their validity beyond leading order of perturbation 
theory observed before. We have thus supported our results for the non-diagonal NLO anomalous 
dimensions of the conformal operators derived in || and we now have a stronger evidence for 
the supersymmetric nature of the universality of the special conformal anomalies of conformal 
operators conjectured there. This issue is currently under study p6| . 

This work was supported by BMBF and the Alexander von Humboldt Foundation (A.B.). 



Appendix. Here we briefly describe the main steps for the construction of an irreducible represen- 
tation of supersymmetry in the basis of the conformal composite operators. The variation of the 
bosonic conformal operators in Eq. (E) leads to 




1 



i-(-i) 





On the other hand since the fermionic operator, (Gip)jj, like the bosonic ones ([|) has to be 
the highest weight vector of the corresponding conformal tower, [{Gip)jj, /C_]_ = 0, with the 
generator of the special conformal transformation /C^, this requirement fixes the indices of the 
Jacobi polynomial to be + — 1. Here d^ and are the canonical scale dimension and spin of 

(2 1) 

the field 0, respectively. For the case at hand we have thus P- . This condition implies severe 
constraints on the coefficients with which the operators (|5|) can enter ^CpOji + G Cj G Oji, which 
form an irreducible representation of the supersymmetry algebra. The corresponding equation for 
*C 'j has two solutions. They are given in Eq. 
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